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Abstract 

We introduce a renormalized IPI vertex part scalar field theory setting in momentum space to 
computing the critical exponents v and rj, at least at two-loop order, for a layered parallel plate 
geometry separated by a distance L, with periodic as well as antiperiodic boundary conditions on 
the plates. We utilize massive and massless fields in order to extract the exponents in independent 
ultraviolet and infrared scaling analysis, respectively, which are required in a complete description 
of the scaling regions for finite size systems. We prove that fixed points and other critical amounts 
either in the ultraviolet or in the infrared regime dependent on the plates boundary condition are 
a general feature of normalization conditions. We introduce a new description of typical crossover 
regimes occurring infinite size systems. Avoiding these crossovers, the three regions of finite size 
scaling present for each of these boundary conditions are shown to be indistinguishable in the results 
of the exponents in periodic and antiperiodic conditions, which coincide with those from the (bulk) 
infinite system. 
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I. INTRODUCTION 



Finite-size effects manifest themselves generically whenever particles or fields are confined 
within a given volume whose limiting surfaces are separated by a certain distance L. Their 
size and shape can affect key properties of the system in comparison with those obtained 
from the L — )■ oo limit ( "bulk system"). Perhaps the most investigated aspects are related 
to critical properties of finite systems l2| , where field-theoretic methods can be employed 
in the vicinity of the phase transitions taking place in the system under consideration. 
Experimentally, the simplest realization of such critical behavior and the role played by 
the finite size corrections show up in parallel plate geometries, for instance, in coexistence 



curves of critical films of certain fiuids 
specific heat amplitudes) in confined ^He 



31 as well as superfiuid transition features (e.g., 
^, |5|. From the theoretical viewpoint, field theory 
studies have been put forth to explain these effects not only for "^He js], but also in thin 
slabs j?, 8| formed by wetting phenomena jsl. The Casimir effect has also been investigated 



finite-lenght cobalt nanowires 



in superfiuid wetting films [10 [.Plus, the recent study of some microscopic properties of 



11| reveals that the infiuence of the finiteness is a ubiquitous 



theme in several properties of physical systems. 

Momentum space e-expansion description 12| of critical properties for finite size systems 



was presented some time ago by Nemirovsky and Freed [NF) [l3|,llj]. The simplest approach 
uses a parallel plate layered geometry, namely, a (slab) volume of material whose limiting 
surfaces (plates) are of infinite extent along {d~ 1) spatial directions and are separated by a 
distance L. This parallelepiped-shaped (e. g., magnetic) material possess a field-theoretical 
description of its critical behavior in momentum space which requires continuous momenta 
components parallel to the (rf — 1) spatial directions and discrete "quasimomenta" along 
the finite size direction of the material. It is basically a combination of effects coming from 
volume (bulk), finite size and surface phenomena. The first two are dominant whenever the 
absolute value of the order parameter (field) is chosen to have the same (not specified a 
priori) value at the limiting plates. (If an external field is allowed in addition to the bulk 
order parameter, and kept at a fixed value at the limiting surfaces, surface effects will become 
proeminent in the discussion of the subsequent criticality, beyond the simpler volume (bulk) 
plus finite size corrections pattern.) 

These geometric restrictions can be realized as many different boundary conditions im- 
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plemented in the bare free propagator. The above mentioned simpler finite size correction 
shall interest us throughout and can be modelled when periodic and antiperiodic boundary 
conditions are employed. (Dirichlet and Neumann boundary conditions mimic free surfaces, 
are appropriate to explain finite size plus surface effects and shall not concern us in what 
follows.) The limitation caused by the boundary conditions provides a scaling variable 7^, 
where ^00 is the (bulk) correlation length of the infinite system. Many computations up to 
first order in e of amplitudes connected to Green functions have been carried out within this 
massive framework, as well as some universal amplitude ratios of certain thermodynamical 

g regions induced by the limitation have 
been proposed. The first one is characterized by 7^ > 1 where perturbative methods can 

Soo 

be applied and the physics is quasi (i-dimensional, characterized by bulk critical exponents 
but limitation dependent amplitudes. The second region corresponds to ~ 1 and it was 
conjectured that the critical behavior is neither (i-dimensional nor {d — l)-dimensional. The 
third region is associated to values of the variable 7^ < 1. It was also argued that in this 
regime the physics is almost {d — l)-dimensional and usual perturbation expansions break 
down jl^. Another prediction stated that the normalization functions and the exponents 
would be the same as those found in the infinite system for the boundary conditions above 
mentioned. 

In this work we introduce a one-particle irreducible (1-P/) renormalized field-theoretic 
version of the NF formalism in order to investigate finite size corrections to normalization 
functions, fixed points, etc., at higher order in a perturbative loop expansion which are 
dependent upon the boundary condition on the plates. Concrete applications for periodic 
(PBC) and antiperiodic (ABC) boundary conditions are explored through the computations 
of the critical exponents t] and u in finite size scaling using the diagrammatic method in 
momentum space, at least up to two-loop order. We improve the understanding of the three 
scaling regions and show that the finite size effects related to the limitation caused by the 
boundary conditions do not show up in the exponents themselves, although they modify the 
ingredients required to compute them. 

We utilize massive fields obeying these boundary conditions on the plates for nonvanishing 
values of L corresponding to fixed finite values of the bulk correlation length. Both first 
region (L > ^00) and the second one associated to finite values of L (— )■ ^00) can be described 
satisfactorily within this massive framework. The remaining region is treated with massless 
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fields having infinite bulk correlation length. In that case, second region is realized through 
the limit L ~ ,^00 — oo. The third region naturally describes arbitrary finite values of 
L and can only be approached using massless fields. The universal results obtained are 
shown to be valid for the three regions determined by the boundedness variable -r- which 
interpolates from infinite to finite (not so small) values of L. The failure of the finite size 
phenomenological scaling arguments regarding the second and especially the third scaling 
region is demonstrated for the first time. 

From our analytical expressions described essentially in terms of elementary primitives, we 
demonstrate that the dominant contribution of the finite size correction goes with the inverse 
power of L only for periodic boundary condition, where dimensional crossover starts to set 
in the critical behavior. This pattern occurs for both massive (t > 0) and massless (t = 0) 
regimes, although with a larger coefficient in the last situation. Antiperiodic boundary 
conditions have the usual crossover at t < in the massive theory as previously discussed by 
Nemirovsky and Freed [l^. Furthermore, our analytical method shows clearly the existence 
of a new type of crossover which takes place for ABC at t > when a term proportional 
to InL becomes important for small values of L. On the other hand, ABC in the massless 
regime t = presents a power law of the type L^^ whenever L is small. Therefore massless 
and massive crossover regimes are completely different for ABC, which is demonstrated here 
for the first time. We show that, as long as we avoid these crossover regions for very small 
values of L, there is no breakdown of the e-expansion into third region and demonstrate the 
validity of the computation of the exponents. As far as critical exponents are concerned, 
the physics of the systems in the three regions is actually quasi d- dimensional, for the bulk 
critical exponents are recovered from the finite size evaluation irrespective of the boundary 
condition and the value of L, i.e., independent of the limitation variable 7^. 

The paper is organized as follows. In Section II we describe the formalism of massive 
fields with fixed finite correlation length (mass). The case L — )■ 00 corresponding to region 
a) is shown to smoothly reproduce the bulk exponents. An introduction to the L — )■ limit 
and how it is related to dimensional crossover is presented as well as the result of the solution 
to higher loop diagrams away from the dimensional crossover region. Section III presents 
the computation of the critical exponents 77 and v using normalization conditions utilizing 
the Feynman diagrams outlined in the previous section, at least up to two-loop level. We 
show that they are L-independent. 
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We set the massless framework in Section IV, using normahzation conditions as well as 
minimal subtraction. We discuss the behavior of the integral for certain values of L and 
compute them in the form suitable for each renormalization scheme. The exponents obtained 
from the setting of massless fields are computed in Section V. For infinite values of L, we 
show that our results correspond to region b). Finite values of L are shown to be equivalent 
to region c). The complete equivalence with the exponents computed using massive fields is 
established. 

In Section VI we discuss our results and point out future potential applications of the 
method to approaching other types of critical behaviors with simple boundary conditions. 
Higher loop Feynman integrals are presented in the appendixes. The massive integrals 
in normalization conditions are described in Appendix A. In Appendix B we display the 
massless integrals in normalization conditions and in minimal subtraction. 



II. MASSIVE FIELDS FOR PBC AND ABC IN THE NF APPROACH 

In this section we begin with a quick review of the NF setup [l^ in order to describe 
our computation of the critical exponents explicitly for periodic and antiperiodic boundary 
conditions. These boundary conditions realize the simplest situation in the discussion of 
finite size effects inasmuch they do not include the effect of free surfaces. 

The layered system can be described by the bare Lagrangian density 

^ = ^1 V0or + ^/io0o + ^Ao0o, (1) 
where 00, yUo and are the bare order parameter, mass (where /Xq = to is the bare reduced 



temperature) and coupling constant, respectively jl6Hl8|. The coordinates are decomposed 
in the form x = (p, z) where p is a (d — l)-dimensional vector characterizing the surface of 
each plate and the z direction corresponds to the region perpendicular to them. The plates 
are parallel and layered in the region between z = and z = L. The field satisfies 0o(^ = 
0) = 00 = L) for periodic boundary conditions, whereas 0o(z = 0) = —(f)o{z = L) for 
antiperiodic boundary conditions. The order parameter can be expanded in Fourier modes 
as 0o(p, z) = Yli j d'^~^kexp{ik.p)uj{z)(j)Qj{k), where k is the momentum vector associated 

j=-oo 

to the {d — l)-dimensional space, Uj{z) are the normalized eigenf unctions of the operator ^ 
whose eigenvalues kj defined by — ^^i^"* = /^"jUjiz) are called the quasi- momentum along the 



oo 

2;-direction. In addition, the eigenfunctions obey the relations ^ Uj{z)uj{z') = 6{z — z') 

j=~oo 

and Jg dzuj{z)u*,{z) = 6jj>. Note that = a{i + r), where a = j = 0, ±1,±2, 
the label r = corresponds to PBC and r = | to ABC. The free bare massive propagator 
(yUQ 7^ 0) in momentum space for either boundary condition is given by the expression 



'Oj y-'J) fc2+<T2(j+r)2+/i2- 

Since a typical Feynman integral involves the product of many bare propagators Gg^- , 
the Feynman rules are modified as follows: beyond the standard tensorial couplings of the 
infinite theory corresponding to a component order parameter, each momentum line 
(propagator) must be multiplied by sj^j^ = dzuj-^{z)uj^{z) and the vertices are multipled 
by the tensor Sj^j^j^j^ = j^^ dzuj-^^{z)uj^{z)uj^{z)uj^{z). Furthermore, each momentum loop 
integral in the finite system can be obtained from the infinite system through the substitution 

oo 

J d'^k — X] ^ / d'^^^k. The eigenfunctions actually depend on r and can be written 

j=-oo 

as u^l\z) = L-hxp{iKjz) which implies S]^]^ = ^h+j2,o and S^llj.^j^ = L-^dj.+j^+j.+j.fl. 
This means that the quasi-momentum is "conserved" along the z direction for periodic and 
antiperiodic boundary conditions. 

Let us define the renormalized IP/ vertex parts from the NF construction. Although 
they do depend on the boundary conditions, we shall not introduce this additional label 
on them. Consequently, considering an arbitrary IP/ divergent (but regularized, say, by 
a cutoff A) bare vertex part including composite operators F^^'^^) {{N,M) ^ (0,2)), the 
statement of multiplicative renormalizability amounts to finding renormalization functions 
zj,^-*, z'^p such that the vertex parts defined by 

^f^\vuHM.9.l^) = (4^))^(4;))^r(^'^^)(p„z„QM;,Ao,/io,A), (2) 

are automatically finite (when the regulator A is taken to infinity). 

In the massive framework, the primitive divergent vertex parts of this field theory are 



chosen to be renormalized in the standard way [16|, but now they are explicitly dependent 
on the boundary condition, even though we omit the label r in all vertex parts. Then, 
we choose the following normalization conditions at zero external momenta (and quasi- 
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momenta), namely 



^Rik = 0, J = 0,g,f^) =fx' + a^T\ (3a) 
grg^(fc,j = 0,^,/i) 

Qj^2 lfc'=o = 1' y^^) 

Tf{h = Q,^l = Q,g,^^)=g, (3c) 

rg'')(A; = 0,j = 0,Q = 0,/ = 0,^7,/i) = 1. (3d) 



Note that the normalization condition on the two-point function above amounts to choosing 
the renormalized mass /i independent of the boundary condition [6|. These conditions are 
sufficient to formulate all vertex parts which can be renormalized multiplicatively. 

First let us discuss the situation at the critical dimension d = A. In that case, utilize 
implicitly a cutoff A to regularization of the integrals and suppose that after the renor- 
malization procedure is defined the limit of infinite cutoff can be taken. We can obtain a 
Callan-Symanzik equation which describes the scaling regime through the following steps: 
i) apply the derivative ^ over the bare vertex part T^^'*^) ((A^, M) ^ (0, 2)) at fixed Aq, A 
which produces the vertex function T^^'^''^^\pi, ii, Qi, i\] 0; Aq, /io, A) at zero inserted momen- 
tum; ii) rewrite the remaining bare vertex parts in terms of the renormalized ones. This 
results in the following expression 

(2p/^ + 4# - - M^)^'r^'"\pi^ Qi, 4 9, /x) = (4) 

where 2p = gg^), ^ = ^^H, f = 4^^, ^ = 4^^. Let the fiow functions 
be defined by the expressions /^^(/i, g){= f ) = ^l'\= = and 7j^(= -p = 



*2 



dlnZ]^- ______ _ _ ,2 



— /i — Qjf~- Multiplying last equation by we obtain its equivalent form in terms of these 
redefinitions as 

(/^l; + - ^lt^ + M^^y)Tf^''\p,^,Q,^',g,^^) = (5) 

2/.^||[4^]-^rif'^^^)(PM.gM;;0,^,/.) . 

Now, taking A^ = 2 at zero external momenta and quasi-momenta and using the normaliza- 
tion conditions Eq.(l3all and Eq. (l3d1l we obtain the Callan-Symanzik equation for finite size 
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systems given by 

+ ^'^^1 - f + ^7i?)rif'"^(P..,Q.4.,^) = (6) 

where 7^^'' = 7^^''[1 + (~)^]- The difference with respect to the original version for infinite 
system is the appearance of 7^""^ which is essentially a scaled version of 7^"^ Note that in the 
limit L —7- 00 ((T —7- 0) we retrieve the original bulk theory naturally within this construction. 
For fixed L, all the discussion of inductive proof of multiplicative renormalizability follows 
exactly as in the bulk system described by the ordinary 0^ theory. 

Away from the critical dimension, the coupling constant has mass dimension. As before, 
let P^'^\g,fi) = be the function which governs the flow of the coupling constant in 
parameter space. In order to get rid of undesirable dimensionful parameters when d = A — e, 
define the Gell-Mann-Low function [f3{g,^)]GL = —^9 + f3{g,f^)- Using the Gell-Mann-Low 
function into the CS equation, it is easy to find that all dimensionful parameters turn into 
dimensionless quantities. For instance, let A = /i^mq be the dimensionful bare coupling 
constant written in terms of the bare dimensionless coupling uo and g the renormalized 
dimensionful counterpart written in terms of the dimensionless renormalized coupling u as 
g = fi^u. Those definitions imply that [l3{g, fi)]GL^ = /^l"")^) i-^-, we get a description 
entirely in terms of the dimensionless renormalized coupling constant, which has a well 



defined scaling limit 19|, |20|. The Callan-Symanzik equation can be expressed in the form 



+ - f + M7j?)rif'^)(p,,.,,Q,,4«,;x) = (7) 

/^^(2 - 70'^^)r5f '■^^^^(p,, ii, Qi, i[; 0, u, /i) , 

where = -e(^), 7^!^) = /^^^H^) and 7^^^ = /^^^H^)- The definition 

zjjp = Z^^ Z^^ can be used to write down another function, namely 7^2'' = P^^^ — g£—) 
which shall be useful to our purposes. The solution of the Callan-Symanzik equation is 
analogous to the infinite systems version and we shall not discuss it here; instead, we shall 
use the results of previous analysis in order to discuss the salient features which naturally 
leads to the ultraviolet fixed points along with the critical exponents for finite systems 
satisfying various boundary conditions. 

Recalling that the infrared divergences are absent in the massive theory, we analyze 
the theory at the ultraviolet region where the momentum of the internal propagators in 
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arbitrary loop graphs are very large, i.e., at the scaling region ^ — )■ cxd [17|, |18|. This means 
that the right hand side can be neglected order by order in perturbation theory just like in 
the field-theoretic description of infinite systems. 

Let us turn now our attention to the computation of the Feynman integrals corresponding 
to one-, two- and three-loop diagrams required to getting the critical exponents r] and u 
perturbatively. 

The one-loop integral contributing to the four-point function is then given by: 

oo ^ 2 



j=-oo 



21 



''[(g + A;)2 + (a)2(j+^ + r)2-f/i2]- 

Remember that fi = = at tree level, where t is the renormalized reduced temperature. 

Performing the transformation p = ^ in all momenta present in the diagram {k' = ^, 

restoring k' — )■ k) and defining r = ^ = (^), we use a Feynman parameter x before 

resolving the integral over p, or in other words, 

oo „i „ 
4^)(fc,z;a,/x) = r/i- dx d'^'^p 

i=-oo>^0 J 

1 

- -2- (9) 

+ 2xkp + xk'^ + r^[(j + r + ix)"^ + x{l — x)P] + 1 



A typical result within our conventions (see Ref. IGj) appropriate to dimensionally regular- 
ized integrals is expressed by the formula 

f d'q _i r{lna-l){m'-k')i-- 

J {q^ + 2k.q + m?Y 2 T{a) ^ ' 

where Sd is the area of the d-dimensional unit sphere. Using this relation we get to 



I',^>ik,t;a,fi) = rfi-^-Sd-iTi^—)Ti2 



2^2'^ 2 



„1 oo 

dx [x{l - x)(fc2 + iV) + r^{j + T + ixf + l]'^"^ (11) 

Jo 



]=-co 



Notice that oc j here is the boundedness variable in the massive theory, where ^ is 
the fixed bulk correlation length. After factoring out the term in the last integral, we 
can proceed in the computation by noticing that the remaining summation can be identified 
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with the generahzed thermal function 21 1 



n=— oo 

= ^f-^ + U<^M (12) 

where 

oo 

/„(a,6) =4^cos(27rma)(^)°-^K„_i(27rm6), (13) 

m=l 

and Ky{x) is the modified Bessel function of the second kind. The identifications a{x) = 



T + ix, b{x) = r ^\/{k'^ + r'^)x{l — a;) + 1 and e = 4 — d permit us to write 

X [r(^)[x(l-x)(P + .V) + l]-t + /.^.(a,6)]], (14) 

Now, using the identity y/7rT(^^)Sd-i = T(^)Sd and expanding in e the argument of the 
gamma function, we can rewrite last integral as 

/f^(A;,i;a,/i) = Sd/i-^[-(l- ^) X / dx[x{l - x){k^ + + 1]--^ 

^ ^ Jo 

+-r-T(2--) / dxfi_^.{T + ix,r-^^x{l - x){k^ + Pr^) + 1)]. (15) 
2 2 Jo ^ ^ 

Whenever we perform a loop integral, the area of the unit sphere 5"^ naturally takes place 
and this angular factor can be neutralized in a redefinition of the coupling constant. We 
adopt this procedure henceforward in all loop integrals and suppress this overall factor. We 
then find 

l![\k,i,r) ^ ^^p^ = ^((1 - f) j; dx[xil - x)ie 

+rH^) + l]-f + |r(2 - ^)Fi^\k,t;r)), (16) 



where 



and 



F^j\k,i-r) = r-^^ dxf.^^{T + xz,h{k,z,r)^, (17) 



h{k, i, r) = r-'^^/x{l- x){k^ + rH^) + 1. (18) 

Since we are going to use normalization conditions in this massive setting, we are interested 
in the simplest situation which occurs for vanishing external momenta and quasi-momenta 
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{k = 0,i = 0). In that case, FY' (r) = r ''fij^t{T,r Furthermore, recalhng that the 
finite size contribution is 0(e°), the one- loop integral can be written as 

I^;\r) ^ lt\k = 0,^ = 0,r) = - |) + ^/.(r,r-i)]. (19) 

In this massive approach, so long as takes finite nonzero values, last equation is suitable 
to compute amplitudes and other observables and demonstrates explicitly that the one-loop 
bubble at the symmetry point is decomposable in the form bulk contribution plus finite size 
correction. 

On the other hand, in the computation of the critical exponents, in practice we use this 
form without any specification of the correction function, since it does not have singular 
behavior for these values of and the divergence structure of this diagram is governed by 
the first term (bulk) in that expression. In that case, the e-expansion is well defined and 
we can proceed to the computation of loop diagrams of primitively divergent vertex parts 
in order to renormalize the theory and obtain those universal quantities. Note that even 
though the last expression is written in terms of /i (t, r^^), we prefer to write the correction 
in terms of F^\k = 0,i = 0;r) anticipating our future discussion of the massless case. In 

2 

addition, applications of the present method might be important to compute amplitudes. 
We would like to understand the importance of the finite size correction in limit values of 
the boundedness variable in calculating an arbitrary amplitude. This amounts to figuring 
out the approach to the regions — )■ (or r — )■ oo) as well as r^^ — )■ oo and what are the 
limits of validity of the e-expansion. 

In order to describe these asymptotic values and its effects on the finite size correction, 
we write the latter in the form 

oo 

F^/\k = 0,i = 0;r) = Ar-''^cos{27TnT){7m)^ K^{27Tnr-^) . (20) 

n=l 

This expression shows clearly that the correction has no poles in e. Since it is well 
behaved, take e = in this whole expression in order to rewrite Eq. ffTOj) at the symmetry 
point as 

lispir) = (1 - |) + 2 J] cos{2ixnT)Ko{27^nr-')\ . (21) 

n=l 

The limit — t- oo corresponds to — t- oo, whereas < < oo represent finite values 
of J. Let us focus our attention in the limit r^^ — ?■ oo. Using the asymptotic form of the 
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Bessel function for x — )■ oo, namely Ka{x) = ^[1 + 0{l/x)], for r — )■ oo one learns 

that the correction term has the behavior 

Urn {^FT{k = 0, i = 0; r)) = ^— V cos(27mr)n-i/2g{-2.nr-i)_ ^^2) 

V / n=l 

Using a trivial unequality to simplify our task, we can show that last term vanishes in the 
wanted limit as follows 

^ oo _ oo 

liyn Vcos(27rnr)n-i/2g(-2.nr.-i)^ ^ ^.^ cos(27rnr)e(-2™^"')) = 

r ^— >oo r- 2 ^ ' — ^— ^-oo r 2 ^ 
n=l n=l 

/^m (- ,\ , . , + ,\ , - 2) ^ 0. (23) 

Therefore, the integral turns out to reproduce the (bulk) value from the massive theory of 



the infinite system [17|. We can identify region j — )■ oo with usual bulk critical behavior 
L — i- 00. The region ^ > 1 interpolates from finite size corrections to the bulk critical 
behavior. 

As J decreases, the finite size correction will increase until it will eventually become as 
big as the pole in e, modifying the leading singularity of the four-point function. To see this 
let us consider the potential trouble which is hidden in the different values of and, in 
particular, in the limit — t- (L — 0). 



22|, 



Let us perform the sum which appears explicitly in the correction term. From Ref. 
the identity 

1 TT °° 1 

1 TT °° 1 1 

"2^^] ^ 2^f^x2 + (2n7r^teP " 2^^^' 

which is valid for positive finite values of the variable x, along with the identifications 
X = 27rr^^, t = TT (and 7 = 0.57721566... is the Euler-Mascheroni constant) implies the 
following result to the one-loop graph 

h{k = 0,^ = 0, r-i) = /X- [1(1 - 1) + 7 + InC-^) + 



T 

00 



+ 2|j[yF^TFT7F n^^ 2^^^^r-^ + (n-ry J^' ^^^^ 
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The simplest way to prove that the two infinite series are convergent in the hmit r ^ — )■ 
is to set directly r^^ = and compute this correction [23]. It becomes 

um (-y \ ] + -y\ ]) = - y\^ - -] + - y\^ - -]. (26) 

n=l n=l n=l n=l 

Now, from the definition of the dilogarithm function ip{l + z) = —7 + Yl'^=i (n+z) ' together 
with the relation ip{l + z) = ip{z) + \ and the value V^d) = —7 — 2Zn2, we easily obtain 

^00 _ 00 

.^^o^sSZM + sEC ]) = (2/n2-l)5,,., (27) 

n=l n=l 

which is finite as advertised. Therefore, for small values of we can write the one-loop 
bubble as 

/2(A; = 0,z = 0,r-i)=/i-^[l(l-^)+7 + /n''' ^ 



■e' 2' ' ' 2 ' 2V7 



-2 + r2 



+ {2ln2 - 1)6,^1 +0{r-^)]. (28) 

This expression for the one-loop four point function depends on the boundary condition. 
It gives the support to identify two types of crossover in finite systems presenting these 
simple boundary conditions away from the critical point (t 7^ 0) as follows. 



First 
Refs. 13 



y, our analytical expression above is a transliteration of the analysis performed in 



14j concerning the breakdown of the expansion in e = 4 — d, namely, when the 



argument of the square root term in the above expression vanishes. Indeed, for periodic 
boundary conditions r = and perturbation theory is invalid in the limit = ^ — > 0. 
For antiperiodic boundary conditions, however, if the temperature is below the bulk critical 
temperature (t < 0), whenever r"^ = —^^^ = 1 the inverse square root blows up. This 
effect was denominated "dimensional crossover" as discussed previously by those authors. 

Secondly, if the value of is decreased further for fixed t > 0, i.e., diminishing L, the 
logarithm term starts to become important for antiperiodic boundary condition when its 
argument becomes around the same order of magnitude that the dimensional pole ^. If 
we switch to cutoff regularization for a moment, the ultraviolet regime is characterized by 
- — )■ ln{—) = ln{A^) with ^ 1. The logarithm contribution will eventually become 

e fj. 

comparable with the ultraviolet dimensional pole, whenever (j) ~ j^^, i.e., when LA ~ 1. 
In terms of a lattice parameter a, A ~ ^ which implies L a. It is the reduction of L for 
fixed t > in this massive framework which is responsible for this new effect. This is a novel 
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type of crossover which only happens for antiperiodic boundary condition at t > and is 
straightforward from our analytical expression given purely in terms of elementary primitive 
functions. This new type of crossover starts when instead of a large number of parallel 
plates, there are only two parallel plates (the limiting surfaces) and the bulk description is 
no longer reliable. Note that this behavior is also there for periodic boundary condition, but 
the square root term proportional to (-1) ^ is overwhelming in that limit. 

A word of caution here. It is dangerous to take the limit t = (or ^ — > oo) in the 
above expression. The reason this limit is inconsistent in this massive framework is that 
the scale invariance of the renormalized theory only takes place in the ultraviolet regime. 
The most appropriate strategy would be to start from scratch with massless fields which 
are scale invariant at this infrared regime, renormalize the theory at nonvanishing external 
momenta scale and push forward all the consequences which follow from this approach. We 
are going to study this case later on and shall prove from a full two-loop calculation that 
the phenomenological scaling theory, which states that e-expansion results have meaningless 
results at i = is incorrect. We postpone this discussion to Sees. IV and V. 

Without loss of generality we can choose /x^ = 1 which is equivalent to a fixed (arbitrary 
but finite) correlation length, such that r = a = (In our subsequent discussion we can 
reconstruct the ^ dependence through its multiplication by L~^.) In fact we could have 
started directly from this choice for the mass scale, and it will define all other massive 
loop integrals yet to be discussed. We have only to keep in mind that this choice makes L 
dimensionless. 

Prom now on, we stay away from the region of crossover in order to compute the higher 
loop integrals. These objects can be computed analogously to our previous one-loop discus- 
sion and the reader is advised to consult Appendix A for details. One typical example is 
the integral contributing to the four-point function at two-loops, namely 



h{ki,k2,h,k4,ii,i2,i3,i4,cr,lJ,) = / -r^-- — ^^rr^- — : ^ 

. ^ J [Qi + cr^iji + 'T 

JlJ2 = -00 

1 



_ y+pi') 

JlJ2 = -00 



X 



(?i + (T'ij2 + t)2 + [(^^ _ 52 + + a^{j, - h + i3 + t)2 + 



""[{P-qir + a^p-ji+rr + l^'] ' ^^^^ 
where P — ki + k2 is the external momenta along the plates and p = ii -|- ^2 is a discrete 
"external" quasi-momentum label. Taking /i — 1 at zero external momenta and quasi- 
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momenta, let the integral be denoted by li^\o, 0; a). From the result computed in Appendix 
A, the outcome to this diagram is 



4k^) = ^((i-^)+e/i(^>0)- (30) 



The integrals contributing to the two-point function at two- and three-loops, respectively, 
turn out to be 

d'^-^qid'^-^q2 



. ,.. + a2(ji + r)2 + ;x2)(g2 + a2(^-2 + ^)2 + ^2) 
[{qi + q2 + kf + + J2 + * + rf + /i^] ' ^^^^ 



and 



J,.. . 2 / d'^-'q.d'^-'q^d'^-'q, 

jlij2 — CO 

1 

(gf + ^x2(j3 + r)2 + [(gi + g2 + A;)^ + + J2 + ^ + rf + /x^] 

[(gi + 93 + A;)2 + a2(ji + J3 + ^ + rY + • (^2) 
In passing, we note that the massless diagrams to be studied later, can be obtained from 
the above expressions by setting /i = 0, although we work in the massive case at = 1. 

The objects required for our purposes are the derivative of those integrals with respect to 
the external momenta, setting the external momenta at the symmetry point in the end of the 
process. Let the derivatives in relation to /c^ be /g''^^ and l'^'^'' (at null k). It is demonstrated 
in Appendix A that they are given by: 

4t(^) = -^(l-| + ^W^^^H^)), (33) 

and 

4%) = -^(l-J + |w^^^H-)), (34) 
where W^^\a) = G^^\a) + H^^\a)-ArFQ^"\a) and F^^^\a), G^^\a) and H^^\a) are defined 
in Appendix A, Eqs. flA7p . (IA14aP and flA14bp . respectively, computed at e = 0. We are not 
going to write them down explicitly since we shall show soon that they will be eliminated 
during the renormalization process. 

With the information at hand, we can proceed to compute the normalization functions, 
Wilson functions, the repulsive fixed point, the anomalous dimension of the field and of the 
composite operator. This task shall be tackled in the next section. 

15 



III. CRITICAL EXPONENTS FROM FINITE SIZE WITH PBC AND ABC 



Now we describe the normalization functions and Wilson functions in terms of the loop 
integrals. These are the fundamental quantities needed to uncover the diagrammatic com- 
putation of universal quantities. 

The occurrence of a nontrivial ultraviolet fixed point, the scaling limit in the ultraviolet 
regime and the simplification achieved for the renormalized vertex parts computed at the 
fixed point are important aspects with consequences in this sort of computation. 

As previously discussed, the ultraviolet flow in momentum space can be described in 
terms of dimensionless coupling constants. In this way, we can write the dimensionless bare 
coupling constant and normalization function Z<^, Z^2 as power series in the dimensionless 
renormalized coupling constant as u^q^ = u{l + o!{\ + a^2^u'^), Z^^^ = 1 + b^2^u'^ + b^^^'u^ and 
Z^2^ = 1 + c[ \ + c^^"^. The divergence structure of these objects are dimensional poles 
appearing as inverse powers of e{= 4 — d) . In order to figure out explicitly each coefficient as 
parameters depending on the loop Feynman integrals computed (at symmetry point defined 
at zero external momenta and quasimomenta) so far, express them in the form 

u- u\l I ^^ + '^ J^-)k I + 

Uo— U[L-\ -'25P^+V 

(iV^ + 6iV + 20)(/gp)^ (5iV + 22)/igp {N + 2)l',P 
1 + + + ^)(^ + ^)(p3?^ - ^).^ (35b) 

(7V^ + 77V + 10)7gp (7V + 2) (7V + 2)(7gpf (,) 

+ [ 6 6 ^-/45pJJM. [60C) 

The flow functions describing the parameter space as the momentum scale varies are 

I3'^^^\u), 'ylp \u) and ^y^J . When they are written as series expansions in terms of u, we 
obtain explicitly 

= -eu[l - a['K + 2{{a[^^f - (36a) 
7^^^ = -eu[2b^;^u + (3^-) - 2b^;^a^;^)u% (36b) 
= 6x.[cS^) + (2cr^ - (cS^))^ - aS^^cS^Vl- (36c) 
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We then employ the results for the integrals presented in the last section for periodic and 
antiperiodic boundary conditions in order to find the values of each coeficcient. 

Looking at Eq. fl36p and comparing it with the expansions of the dimensionless bare cou- 
pling contant, normalization constant of the field and that from the composite operator, we 
can read off their results. The fact of the matter is that the utilization of normalization 
conditions provokes the appearance of correction terms in those functions which depend 
explicitly on the boundary conditions. As is well known, this is a prevalent artifact taking 
place in this renormalization scheme. Thus, the renormalization functions at two loops will 
not be equal to those provenient from the bulk system. Nevertheless, these nonuniversal 
corrections are going to cancel out in the expression for universal quantities as we shall see 
in the remainder of this section. 

Working out the details using this prescription, it is not difficult to prove that they are 
given by the following expressions 



a 



+ 1 1 



(r) 



[l-7Te+7T/i(r,a-i)e], (37a) 



=(-6^) 36^ + ( 36^ ^ 

+ + (37b) 

hf = -^^[l-\^m^^^a)l (37c) 

4^^ - Sl^mip^y^ _ !! + (.,.-)], (37d) 

^ 1296e2 L 4 2''2^ ' ^J' ^ ' 

^ (iV + 2)(^ + 5) _ + 17iV + 26 (iV' + 7iV + 10) 

2 36e2 72e 36e ''^^ ' ' ^ ' 

Using the coefficients well as into Eq. fl36ap we get to 

= + ^-^(1 - 2 + 2/|(^'^ ))« ^« • (38) 

Furthermore, we can also find the solutions for the functions related to the anomalous 
dimension of the field and composite operator. They read 

if = - i + - ^^(1 + iy'"(<^) - /i(r.<^-'))«l. (39a) 

,(?.(i^^(l-l + l/.,.,.-))„-(^„^ ,39b) 
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The eigenvalue condition /3^'^^{u^) — yields the repulsive ultraviolet fixed point which is 
given by 

As usual we identify the fixed point value of the anomalous dimension of the field with the 
critical exponent rj, i.e. 

Similarly, it is straightforward to show that 

(^), , (iV + 2)e.^ (6iV+ 18) , 

which in conjuminance with the equation ^^(uoo) — 2 — rj — yields the value of the 
correlation lenght exponent 

1 (iV + 2) (iV + 2)(iV^ + 23iV + 60) , 

2 4(Ar + 8) 8(iV + 8)3 ^ ^ 

Note that these exponents are indeed independent of the boundary conditions and are ex- 
actly the same as those obtained from the bulk system through an analogous utilization 
of diagrammatic methods. If we use the remaining bulk scaling relations, we find that all 
the critical exponents for these simple boundary conditions reproduce those from the bulk 
confirming the one loop analysis by Nemirovsky and Freed and extending it to the present 
higher loop correction. 

We have thus succeeded in formulating IPI renormalized vertex parts for the massive 
theory in order to compute the exponents within the perturbation expansion in e. We shall 
introduce the framework of massless fields in the calculation of the critical exponents in the 
next section. 



IV. NF METHOD FOR PBC AND ABC USING MASSLESS FIELDS 

We start by describing the method of renormalization for massless fields, given by the 
previous Lagrangian density. The (d — 1) -dimensional momentum space lies along the plates 
and characterizes directions perpendicular to the finite size direction L, which is represented 
by the quasimomenta. The free bare critical (massless) propagator is G^j(k, j) — ■ 
The former definitions for the tensors and Sj,i^j„j. hold for the massless case as well. 
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Let us focus now in the renormalization scheme to be chosen in the massless situation. 
The massless theory has infrared divergences at zero external momenta, so it would be wor- 
thy to define the renormalized theory of the few bare primitively divergent vertex parts to 
be finite after using the divergent normalization constants. It is very simple to employ two 
independent renormalization schemes in order to compute critical exponents: normalization 
conditions and minimal subtraction. We emphasize that the integrals involved can be re- 
solved for arbitrary external momenta so that the most convenient form of the results can be 
used to pursue the computation of the critical exponents in either renormalization scheme. 
We postpone this discussion to the next section. For the time being we shall analyze ba- 
sic facts regarding the structure of the normalization conditions. In the present section, a 
thorough analysis of the one-loop four-point diagram in different limits shall be worked out 
along with the quotation of higher loop diagrams (extracted from Appendix B). 

Normalization conditions is appealing due to its simplicity after the choice of the sym- 
metry(subtraction) point for the external momenta taken at nonzero value. Let ki be the 
external momenta of a (d — 1) -dimensional transversal space and let k be the external mo- 
mentum scale along the plates where the renormalized theory is defined. At the symmetry 
point we choose ki.kj = ^(45jj — 1) leading to {ki + kj^ = k^. We fix the external momen- 
tum scale of the two-point function at /c^ = ac^ = 1. The multiplicative renormalization can 
be achieved through conditions on the primitively divergent bare vertices at zero mass, such 
that their renormalized versions take the following values at the symmetry point: 



It is important to mention that the symmetry point implies that the insertion momentum 
in last equation satisfies — {ki + k2y. 

Multiplicative renormalization arguments can be most easily implemented when the bare 
theory is regularized through the ultraviolet cutoff. Indeed, when the normalization condi- 
tions given above are replaced into the renormalized vertex parts defined by 




(44c) 
(44d) 



(44b) 



(44a) 



r^'^'^^iv i O 

R yynt '"n-i ^;; 



4,^,0) = (z(:))f(z(^))^r(^'^)(p„,i„,g. 



i^, Ao,0, A), 



(45) 
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they turn out to render them automatically finite when the regulator A is taken to infinity. 
Although we mentioned the cutoff, we could also use another regularization scheme. In fact, 
we shall shift the argument to consider dimensionally regularized diagrams. We shall utilize 
the cutoff whenever we can successfully simplify the point under consideration. 

Imposing that the bare theory does not depend on the momentum scale where the renor- 
malized theory is defined we find a renormalization group equation in terms of dimensionful 
quantities. At the critical dimension the coupling constant is dimensionful just as discussed 
in the massive setting. Away from the critical dimension, similar arguments can be devised 
to go from dimensionful quantities to dimensionless amounts. 



Let the flow functions be defined by the expressions f3^'^\K,,g) = n^, 'j^J^ = k— ^ 

and = —f^ — 'Bn' renormalized (bare) dimensionful coupling constant is defined in 
terms oi k g — k^u (A = k^uq), where u {uq) is the dimensionless renormalized coupling 
constant. In order to get rid of undesirable dimensionful parameters when d = 4 — e, 
define the Gell-Mann-Low function [^{g, k)\gl — —eg + ^{g, i^)- Consequently, we find that 
[(5{g,K)\GL§^ = /^('w)^) ^-iid the renormalization group equation for the multiplicatively 
renormalized vertex parts read 

{>^Yn + - f ^i^^ + M^V)^f^\Pn, in, Qm, C 0) = 0, (46) 

where ^W(«) = -e(^), ^,'^;\u) = and 7^^ = /S^^H^)- The combina- 

tions Z^2^ = Z^2Zj^^ and 7^2^ = /3'''^\ — shall also be employed. We emphasize that 
the dynamic variable now is the external momentum scale where the renormalized theory is 
defined. The solution is identical to that from the ordinary 0^ theory and it is not going to 
be discussed here. As discussed in the massive theory, we employ solely the definitions above 
for the sake of determination of fixed points and other universal quantities via diagrammatic 
tools. 

To begin with, we write down the one-loop contribution for the four-point function, 
namely 

1 



""W+W+WU + i + W]' ^^^^ 
We utilize Feynman parameters to solve the integral over the continous momentum using 
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standard formulae for dimensional regularization. Using Eqs. (fTOj) and (fT2i) along with the 
identity y/7iT{^^)Sd-i = T{^)Sd, where Sd is the area of the ci- dimensional unit sphere, we 
obtain a result proportional to Sd- As before, we neutralize this angular factor appearing in 
each loop integral by absorbing it in a redefinition of the coupling constant. Collecting this 
steps together we determine the massless expression for the four-point one-loop contribution 
in the form 

/^^ (A;, a) ^ = i ((1 - f ) Jo dx[xil - 

+aH^)-^2 + |r(2 - l)Fi^\k,i;a)), (48) 



where 



and 



Fi^\k,i-a) = a-^'' dxfi^^(T + xz,h'{k,i,a)^, (49) 



h'{k, i,a) = a-^ ^/x{l-x){k^ + aH'^). (50) 

Note that Eqs. (H^ and (l5Up are the massless counterparts of the massive definitions Eqs. 
ffT7|) and f[T5]) . respectively. Using the representation flT^ . the above definitions lead to 

Fy\k,i;a) = Aa-'Yl Jq dxcos{2TTn{T + ix)){- 



n=l [^(l-^KfcW^ 



xir|(27rn(T-i[x(l -x)(P + (51) 

For both normalization conditions and minimal subtraction, the external quasi-momentum 
label can be taken as the zero mode value {i = 0) without loss of generality, which simplifies 
our task. Recalling that the finite size correction is 0(e°) and neglecting 0(e) terms, we can 
rewrite Eq. (H5]) as 

I^^\k,i = 0;L) = i(l - f - f dxln[x{l - x)^]) 

°° 1 1 
+2E Jo dxcos{27mT)Ko{nL[x{l - x)k'^]2). (52) 



n=l 



We can readily take the limit L — )■ oo (a = -^) by considering the Bessel function for 
large values of its argument. Its asymptotic value is given by Ka{z') = y^^e^^' [1 + O {1 / z')] 



with z' = nLk^Jx{l — x) = nB and n = 1, 2, .... Therefore, the correction becomes 

/ dxcos{2nnr)KQ{nL[x{l - x)k^]^) = \ — dx[x{l - x)]~^ 
n=i Jo V Lk Jo 

oo 

'^^cos{27mT)exp{—nB). (53) 



oo 

X 

n=l 
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The remaining series is a geometric one that can be computed when we express the cosine 
in terms of the complex exponents. We can write it in terms of an upper bound through 
the inequahty 



oo ^ 

^cos{2.nr)exp{-nB) = _ ,,^(_^)^^^^ ^ ,,^(_^)^^^^ " l] < 1- 



(54) 



If we take directly the limit L — t- oo before performing the integral we see that this term 
vanishes and multiplies the prefactor which is also zero. Instead, if we use the upper bound 
we can estimate the integral which in the limit L — oo implies the result 

^r[^(|)]^ 



oo I 

2[E !ldxcos{2'KnT)Ko{nL[x{l-x)e]-^)\ < 



-> 0. (55) 

Therefore, the finite size correction interpolates from the contribution for large but finite 
values of L and vanishes to infinite values of L even in the massless case. Sometimes, it is 
useful to define the bulk parametric integral = dxln[x{l — x)/c^], take k = \fW and 
consider the correction term at e = 0. 

Let us try to comprehend the limit L — )■ 0. Using Eq. flMj) with the replacements x — )■ 



z = La/x(1 — x)A;2, t = ^ and taking into account the above observations, we find at i = 
a simple expression useful for minimal subtraction 

li-\k, . = 0; L) = 1 (1 - I - U{k)) + b + rfx/n[^^^Mlzf))i]] 
If' 1 1^ , 1 1. 



dx — H — > dx 



, oo „i 



(t' + x{l - 2,„i-'« \/((n-r)2+i(l-i)||i]2) " 

OO „i 

K=iJo ^((n + r)2 + x(l-x)[f]2) ^ 

whereas computing it at the symmetry point convenient for normalization conditions (fc^ = 
1), we obtain 

/£),(.) = 7 (1 + |) + [7 + 

1 1 r' 1 1 

+ - / dx[ — ] 

^(r2 + x(l-x)[^]2) K=iJo ^((n-r)2 + x(l-x)[^]2) ^ 

^ OO 1 

+ 2 E / / - -]• (57) 

K=Jo ^((n + r)2 + x(l-x)[^]2) ^ 
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Except for the integrals, the L-dependence in last equation has pretty much the same 
form as Eq. fl25l) has in = ^ described in the massive case. In that case the theory 
renormalized at "mass" ^ = 1 (fixed and finite bulk correlation length) and zero external 
momenta is completely analogous to our renormalized massless theory at the symmetry 

= 1 (infinite bulk correlation length). After those choices we just have to recall the 
variable L is dimensionless. Except for minor modifications like the extra integral on the 
Feynman parameter x the discussion of the various terms parallels that for the massive case. 

In order to study the limit L — )■ 0, let us start with the last two integrals. It is licit to 
take L = inside both of them, such that their summation produces precisely Eq.( l271) . The 
first integral is straightforward and its result is /^[^] — 1- We are then left with the task 



of evaluating the second integral. The identity 22 1 



dx —1 r2cx + b-. 

arcszn[ — J , (58) 



Va + bx + cx^ V-A 
is valid for c < and A < 0, where A = 4ac — 6^. Performing the identifications a = and 
^ = = set to 

dx 47r . r 1 -, , 

= = —arcsin[ — =J . (59) 



Collecting these steps together, we can rewrite the one-loop contribution to the four-point 
function in the form 

= ;(1 + ^) + b + (2^^2 - 1)5,,. + ln[^] - 1] 

27r r 1 -| , . 

H arcsinl — = . (60) 

Last equation shows that last term inside the finite size contribution for periodic boundary 

2 

condition r = becomes exactly ^ which differs from the expression in the massive case 
(with fixed ^ = 1) by a factor of vr due to the effect of performing the integral over the 
Feynman parameter x. The new situation occurs for antiperiodic boundary condition (t = |) 
in the last term, whose limit L — becomes -p-. The logarithm is still there just as before in 
the massive case, but now these two different power-law behavior present in both boundary 
conditions are dominant in the L — limit. Therefore, the finite size correction is generically 
enhanced for both boundary conditions in the critical massless theory consistent with the 
enhancement of fiuctuations at this regime. At the critical point crossover starts earlier 
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in ant iper iodic than in periodic boundary condition as evidenced by our analytical result 
in above equation. Thus, the massless case has nontrivial aspects in comparison with the 
massive case, as shown here for the first time, since the second type of crossover discussed 
in Sec. II for t > in antiperiodic boundary condition is now absent for t = 0. 

Provided we stay away from the crossover regions characterized by very small values of 
L, the two descriptions are almost equivalent, even though the equivalence is not complete, 
as far as the finite size contribution is concerned. The e-expansion is well defined in both 
situations, if the crossover regions are precluded from our analysis. What is really remarkable 
from the massless and massive analysis is the bulk correlation length independence of the 
finite size correction. From last equation, finite values of L persist in the correction even 
when the starting point of the massless theory corresponds to -I — )■ 0. Thus, region c) is 
available to our scrutiny. Consequently, the previous phenomenological conjecture that the 
massless limit cannot be understood in terms of e-expansion is unfounded. From now on, 
we are going to consider finite (but not too small values) of L, since the finite size correction 
has a good behavior and e-expansion methods can be utilized without further problems. 

We leave the task of computing the higher loop integrals to Appendix B. Here we simply 
quote the results. The massless counterpart of the integral which contributes to the four- 
point function at two-loops can be extracted from Eq.f l29l) by setting /i = 0, namely 

1 



{ql + ^'(J2 + tY) [(gi -q2 + + a\3, - 32 + ^3 + rf] 

''[{P-q.Y + a^p-n^Tf] ' ^^^^ 

At zero external quasimomenta, there is no loss of generality to approaching the minimal 

subtraction scheme with arbitrary external momenta or normalization conditions with fixed 

nonvanishing external momenta. From our discussion in Appendix B we obtain 

1 , pj 

I^/\kAcr) = ^(1 - 2 - ^^(^) + ^^i^H^'O))- (62) 

At the symmetry point = = 1 it takes the simpler form by setting e = into the finite 
size correction 

4U^) = ^(l + | + e^o^^^(^))- (63) 
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Analogously, the integrals contributing to the two-point function at two- and three-loops 
can be read off from Eqs. (1311) and fl32|) at /i = 0, respectively, see Appendix B. I3{k, z = 0; cr) 
in a form appropriate to minimal subtraction reads 

(64) 

In normalization conditions the derivative of this integral with respect to k"^ computed at 
the symmetry point k"^ = 1 can be written as 

4t(a) = -l(l + |-2eiyo(a)). (65) 

Using a similar reasoning, the solution for J5 appropriate within minimal subtraction of 
dimensional poles was found to be 

lt\k, a) = -l^{{eWT')[l+^-3ez,{k'+a\')]-3eFi:\k, z = 0; a)-6eF|,i(A:, t = 0; a)) . 
be z 2 ^ 

(66) 

On the other hand, its derivative in relation to k'^ at the symmetry point A;^ = 1 is required 
when applying normalization conditions as our renormalization scheme. Hence, 

4p(^) = (1 + 2e - SeWoia)) . (67) 

Note that we can safely replace the values of the subscript of the additional functions ap- 
pearing as finite size corrections at e = into Eqs. f l62p . f l64p and (IB^ when we employ 
the minimal subtraction scheme. This will facilitate the computations of the normaliza- 
tion functions, since many different functions become identical at e = 0. We have now all 
the required integrals to compute critical exponents in normalization conditions or minimal 
subtraction, to which we turn our attention next. 

V. CRITICAL EXPONENTS FROM THE MASSLESS APPROACH 

The results displayed in last section can be substantiated in the calculation of the crit- 
ical exponents u and 1] at two- and three-loop order, respectively, through diagrammatic 
(perturbative) methods. We are going to compute the critical indices and show that they 
are independent of the renormalization scheme, using either normalization conditions or 
minimal subtraction. 
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A. Normalization conditions 



As before, we write the dimensionless bare coupling constant and normalization function 
Z^, Z^i as power series in the dimensionless renormalized coupling constant as Uq^^ = u(1 + 



= 1 + foW^^ + and Z^J = 1 + cS'^M ■ 



c^^"^. A considerable labor 



can be saved by noting that the structure of these equations are the same as those previously 
discussed in Sec. Ill for the massive case. Indeed, using Eqs. fl35|) from Sec. Ill and replacing 
the values of the massless integrals computed at the symmetry point, we can determine the 
above coefficients in the massless theory as well. Following this trend, we encounter the 
following values for the coefficients 



02 



.AT + 8,2 2A^2^23A^ + 86 ,N + 8.2F^ 
' + + i—^) - 



6e 

(iV + 2) 56 

144e ^ 4 
{N + 2){N + 8) 

1296e2 



72e 



2eWt\a)], 



5e 



6e ^ 2 2°^^^' 
(A^ + 2)(A^ + 5) 2N^ + 11N + U (A^^ _^ 7A^ + 10) 



36e2 



+ 



72e 



+ 



36e 



(68a) 

(68b) 
(68c) 
(68d) 
(68e) 
(68f) 



Now, introducing the coefficients a^^^ as well as a^2^ appropriate to the massless formula- 
tion into Eqs. fl36l) from Sec. Ill, we first determine the flow function in the zero mass limit 
which yields 

P^^\u) = -en + (^)(1 + i + '^Ft\a)y - ^-^^u\ (69) 
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The attractive nontrivial infrared flxed point at two-loop order is found out from the condi- 



tion /?( 



u 



0, namely 



6 



u 



9A^ + 42 



1 
2 



)>]■ 



(70) 



Let us now utilize Eqs. fl39ap and (]39bp in order to determine the Wilson functions for 
the massless formalism. Inserting the coefficients determined above into that equation and 
proceeding along the same lines as before, we get the expressions 



+ V - - ^^(1 - 4^o^^V) - 2Fr(a))n], (71a) 



14> 

(.) _ (iV + 2) 
V ~ 6 ^ 



4 

- + - F^^) 
2 + 2^° 



12 



(71b) 
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Computing these functions at the nontrivial infrared fixed point u*, we find that, i) rj = 
'jlj^\u*) is identical to Eq. fHT]) and ii) 'jjpiu*) is equal to the expression in Eq. fH2|) . which 
leads to the exponent u from Eq. fl43p . 

The massive treatment in the ultraviolet regime is therefore completely equivalent to the 
massless framework at the infrared region, for they originate the same critical indices, even 
though the intermediate steps are completely distincts in the two formalisms, in compliance 
with universality. This thorough treatment of massless fields shall be concluded in a moment 
with the computation using minimal subtraction. 

B. Minimal subtraction 

Here we are not going to calculate explicitly the critical exponents. Instead, we are going 
to calculate the fixed point as well as the functions 7^^^ and 7^2'' at the fixed point. As these 
functions at the fixed point are universal, they should be equal to the ones obtained using 
normalization conditions, leading to the same exponents in either renormalization scheme. 

The dimensionless bare couplings and the renormalization functions are defined in mini- 
mal subtraction by 

00 

ul^^=u[l + J24\e)u% (72a) 

i=l 

00 

i=l 

00 

Z^.^ = l + J2^'^\e)u\ (72c) 

i=l 

The renormalized vertices 

T^^\k,u,^) = zj^^T^'\k,u^;\^), (73a) 
T^i\k,,u,K) = (4^yr(^)(A:„4^\K), (73b) 
Tf'\k,, k2,p; u, k) = 4?r(2'i)(A;i, A;^,^, k), (73c) 

should be finite when e — )■ to any desired order in u. Observe that the external momenta 
into the bare vertices are mutiplied by and all the external quasimomenta of the diagrams 
are set to zero in order to simplify matters. The coefficients a['^\e), h^l\t) and cf\e) are 
obtained by requiring that the poles in e be minimally subtracted. The bare vertices can 
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now be expressed as 







Air) (r) 
^2 ^0 



(A 



[u 







(r) (1) 
2 



.(^) (1) 

-'2 



'2 



.(r) (2)w„ (r)x2 



(74a) 
(74b) 
(74c) 



We recognize that 3^"* is proportional to the integral I^' and B'^' is proportional to I'^' . In 
the remainder we shall suppress the upper indices in the integrals referring to the boundary 
condition, but keeping them implicitly. Explicitly, the coefficients can be written in terms 
of integrals like 
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B. 
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(2) ^ (SAT + 22) 

54 

_ (Ar + 2) fc, 



2/^2 + h 
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(2) 



108 
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(75a) 
(75b) 
(75c) 
(75d) 
(75e) 
(75f) 

(r5g) 

(75h) 



Firstly substitute Eqs. ([75]) inside Eqs.(l74]). Next, utilize Eq.fl72al) into Eqs.(l74]). Finally, 
impose that the renormalized vertex parts expressed as Eqs. (1731) are finite via minimal 
subtraction of dimensional poles. Interestingly, all the logarithmic integrals in the external 
momenta as well as the finite size corrections appearing in I2, 13, and J5 cancell each other 
in the algorithm of renormalization. This results in the determination of the coefficients in 
minimal subtraction, or in other words 



U 



7W 



1 - 
1 + 



6e 

N + 2 2 



144e 
N + 2 

6e 



u' + [- 



36e2 
{N + 2){N + 

1296e2 



24e 

(iV + 2)(iV + 8) 3 
+ 5184e J ' 



^ (iV + 2)(iV + 5) (iV + 2) 2, 

U + [ r^-^^ — \U 



36e2 



24e 



(76a) 
(76b) 
(76c) 
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What is amazing is that in minimal subtraction the above mentioned renormahzation func- 
tions do not depend on the boundary condition explicitly, since all that dependence cancelled 
out naturally, i.e., they do not appear in the right hand side of the above equations. 
Furthermore, from the renormahzation functions one can obtain: 

(.) _ iN + 2) _ (iV + 2)(iV + 8) 
7^ - u ^^28 " ' ^^^^ 

(iV + 2) 1 

%2 = u[l--u]. (78) 

The fixed point is defined by (3^^\u*) = 0. Then, we find: 



(9A^ + 42) 
(8 + iV)2 



(79) 



Substitution of this result into the renormahzation constants will give at the fixed point 
, ''^^ = ?7, whereas, in addition, we have 

(.) _ (iV + 2) 6(iV + 3) 



This expression is equal to Eq. fl42|) obtained using the massive method and consequently 
lead to the same exponent u from Eq. (H3l) . Notice that Eqs. fl76l) - (!80l) are the same as their 
counterpart obtained in minimal subtraction for the usual bulk theory. All the construction 
of renormahzation schemes developed in the present work are thus consistent with univer- 
sahty, which states that critical exponents (among other universal quantities) are scheme 
independent as we have shown herein. 



VI. DISCUSSION OF THE RESULTS AND CONCLUSION 

We have computed critical exponents at higher loop order from finite size layered systems 
subject to periodic and antiperiodic boundary conditions on the hmiting surfaces of the slab 
(parallel plate) geometry by defining a one-particle irreducible (1-P/) vertex parts formalism 
to the previous field-theoretic framework of Green functions for those systems introduced by 
Nemirovsky and Freed earlier. In order to do that, we determine normalization functions as 
well as fixed points and show that they depend on the boundary conditions whenever we use 
normalization conditions either in the massive or massless methods. In minimal subtraction, 
however, we find that these quantities are independent of the boundary condition. 
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We confirm that for large values of L and in the L — oo all the finite size corrections 
are under control and the critical exponents obtained in this way are identical to those 
from (i- dimensional (bulk) universality class. In the case of periodic boundary conditions 
we proved that dimensional crossover only occurs at very small values of L, the behavior 
of the finite size correction is proportional to in the massive and massless cases, being 
independent of the value of the bulk correlation length, although the coefficient of thies term 
is larger in the massless case. This extends the previous analysis by NF performed solely 
in the massive case. As far as the crossover regimes are concerned, antiperiodic boundary 
conditions do not present the simple behavior from PBC. In addition to the "dimensional 
crossover" previously discussed by NF below the critical temperature t < 0, we have found 
a new regime of crossover for ABC which exists only for t > characterized by much 
shorter values of L than its dimensional crossover counterpart occurring in PBC, i.e., when 
the lattice constant is of the same order of magnitude of the distance between the limiting 
plates. Moreover, the crossover in the massless case t = for ABC starts earlier than in the 
massive case. Actually, the finite size correction is proportional to in that case, which 
starts to modify the bulk critical behavior for larger values of L than its counterpart in the 
PBC case. Thus, fluctuations at the critical point enhance the effect of crossover in finite 
systems in these simple boundary conditions. 

Let us discuss the connection between one previous two-loop calculation using NF for- 
mulation with our work of higher loop integrals computed in whole detail in Appendix B. 
Actually, Krech and Dietrich Ref. j3] used NF formulation for massless fields in their com- 
putation of Casimir amplitudes (see Appendix A therein). However, only one two loop 
diagram was computed for the free energy, which actually corresponds to a squared tadpole 
diagram. But this is equivalent to a one loop computation, since normalization constants, 
fixed points, etc., were computed only at one- loop order. Unfortunately, perhaps the lack 
of a better representation for the function fa{a,b) at the time of the writing of that paper 
prevented them to obtain simple answers in terms of elementary primitive functions. Con- 
sequently, they abandoned NF method in momentum space in their subsequent work with 
their collaborators. The work presented here, on the other hand, permits us to go beyond 
the simple conclusions of previous analysis: even though the exponents are identical to those 
from the bulk system when we avoid the crossover region which is certainly not too exciting, 
the crossover regimes assessed by the analytical results described in the present paper sheds 
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new light on the fundamental difference between massive and massless regimes to finite size 
systems criticality, and how fluctuations enhance the effect of finiteness in the latter. More- 
over, a consistent description in terms of massless fields with inequivalent crossover regions 
in comparison with the massive case for ABC is certainly a step forward which cannot be 
underestimated. Since the consistency of the critical regime for both boundary conditions 
implies that phenomenological scaling relying on the failure of e-expansion results in the re- 
gion ^ — 7- is incorrect, the present study should be considered the starting point to widen 
the subject and put it on new grounds, such as computing higher order universal quantities 
like amplitudes in order to improve our present knowledge of critical finite systems. 

Away from the crossover regimes, we have shown a complete equivalence between the 
formulations using either massive or massless fields, where the renormalized mass scale plays 
the analogue role of the external momenta scale used to fix the (new) massless theory in 
normalization conditions. A step further is the minimal subtraction treatment for massless 
fields. 

Thus, our critical exponents results confirm the previous expectations pointed out by 
Nemirovsky and Freed concerning a behavior identical to those describing the bulk system in 
higher order loop computations. Contrary to previous conjectures, it is not the boundedness 
variable ^ < 1 which makes the e-expansion invalid, but small values of L decreasing below 
a given threshold which are responsible for crossover. This crossover description is far 
from being completely figured out, but the resources developed in our trend here should be 
encouraging to tackle this problem. We hope our discussion in the present work can serve 
as an introduction to this subject and might be valuable when new perturbative methods 
to treat the limit L — )■ become available. 

A rather interesting topic is to extend the IP/ method at two-loop order to treat systems 
with Neumann and Dirichlet boundary conditions, since they are more appealing from the 
phenomenological viewpoint. They characterize free surfaces 2J|, which disturb further the 
system due to the breaking of translational invariance along the finite directions. First, 
the quasi-momentum are not conserved in the treatment of these boundary conditions. 
Consequently, in order to renormalize the theory we have to introduce distinct external 
fields, one in the bulk and other in the limiting surfaces. This implies that the surface 
parameter becomes relevant and requires a new normalization function to renormalize it. 
In spite of these additional aspects, we expect that this topic can be investigated along a 



31 



similar line of reasoning to that employed in the present work. 

Another intriguing; perspective is to consider the finite size approach to competing systems 



of the Lifshitz type 25l-l27j. It remains to be seen if the competing axes with arbitrary 



momentum powers permit exact results when the finite size direction points along any of 
them. The last few years have witnessed promising new applications of this kind of field 
theory from quantum field theory to quantum gravity and cosmology. The most direct 
application is to study aspects of space(time)s with one compact spatial dimension, called 
"Lifshitz space(time)", e. g ., in the Horava-Lifshitz theory of gravity 28|| and other simpler 



quantum field theories 



29|. 



Finally, it is possible that the results obtained in the present paper must be used to 
update certain computations of amplitude ratios of certain thermodynamical potentials. 
Other aspects like extension of the present method in the analysis of semi-infinite systems 
are also worthwhile. 
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Appendix A: Higher order massive integrals in dimensional regularization 

Since the relevant one-loop integral integral contributing to the four-point was discussed 
in detail in the main text, we shall discuss only two- and three- loop diagrams, using the 
one-loop result extracted from the text whenever possible. 

The required integrals are computed at unity mass, which makes r = ^ = ^ . The 
simpler contributions come from the two- and three-loop diagrams of the two-point function, 
respectively, given by the following expressions 

d'^-\id'^-^q2 



l! + cT^Ji + r)2 + 1) (gi + a^U2 + rf + 1) 



^ (Al) 



[(gi + 92 + kf + +J2 + 1 + r)2 + 1] ' 
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Ji j2=-oo 



1 



(g| + ct\j3 + r)2 + 1) [(gi + g2 + A;)^ + cr^{j, + j 2 + ^ + r)^ + 1] 

[(gi + + + + J3 + ^ + r)2 + 1] • ^^^^ 
The remaining integral is the nontrivial contribution to the four-point function at two 
loops, namely 



Jl j2 = -00 

1 



(g? + ^'(ji + r)2 + i) 



X 



{ql + ^'(J2 + ry + 1) [(gi - + k^f + a^O^ - + 23 + r)^ + 1] 



''[{P-q.f + a^p-n + rY + l] ' ^^^^ 

where P = fci + A;2 is the external momenta along the plates and p = ?i + ^2 is a discrete 
"external" quasi-momentum label. For convenience, we shall compute all the integrals with 
all external quasi-momenta set to zero. 

The systematics to solve J3 and is very similar: in the former there appears a four- 
point one-loop subdiagram, whilst in the latter a squared of that object takes place. We 
first solve the internal bubble(s) belonging to I3 (h), and use Feynman parameters to solve 
the external resulting bubble. It is important to divide each loop integral by the unit area 
of the (i- dimensional sphere 5*^, a standard procedure within this technique. 

The objects required to our purposes are the derivative of those integrals with respect to 
the external momenta, setting the external momenta at the symmetry point in the end of 
the process. 

Let us apply this general strategy to perform a detailed computation of l'^ at zero external 
momentum and quasi-momentum. First, set the external quasi-momentum i = inside 
Eq.( 1Alll in order to simplify our task. Let us rewrite this integral in terms of the one- loop 
subdiagram l2'^\k,i = 0,r) in the form 

/r(M ^ 0;.) ^ . f j,-.,^^_jg-LM^^ (A4) 

Sometimes, it is appropriate to perform the change of variables q' = q + k, such that the 
external momenta is exchanged to the external subdiagram; see below. Now using Eq. (|T6]) 
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along with the definitions Eqs. (fT7|) . ( !T8|) given in the text in order to solve the internal 
bubble, we find 



it^ (k, ^ = 0■,a) = al ((1 - f ) J^' dx Ef^^ / 



+ f r(2 - t) / + k,^ = 0■, a)) . (A5) 

Before going ahead define the objects 

F'(^)(.) = ?^2(^| (A7) 

i3(A;,a,x)= V / — r, (A8) 

which shall be important in what follows. Inserting last equations into the expression for 

It) 

, the latter can be rewritten as 

4''\k,i = 0;a) = dx[x{l - x)]~U3{k,a,x) 

+s/-r{2-'-)F^]l{k,t = 0;a)y (A9) 

Next, let us compute is^k, a, x). First, perform the change of variables q' = q + k. We utilize 
another Feynman parameter which leads to 

X ^ . (AlO) 



- 2ykq' + yk^ + y + a'^{j + yrf + y{l - y)aH'^ + -^^^^ 



Resolving the momentum integral using Eq. ljTOi) . we obtain after expanding m. e = A — d the 
result 

g,_ir(^)r(-i + e) f\ 

t3{k,a,x) = ^ / dy{l-y)^ 



2T{ 



2> -JO 



1 — — e 

^ [^(l_^)[P + ^V]+a2(j+2/r)2 + 2/ + ^^^:^]'" . (All) 



j=-oo 
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Utilize the identity ^/7lT{^-^)Sd-l = T{^)Sd, absorb 5*^ in the redefinition of the couphng 
constant just as explained for the one-loop four-point graph in the main text and expand in 
e. Pluging these steps in I^^\k,i = 0; a), it is easy to show that 



e r(2-f)r(-l + e) 



dx[x{l — x)] 2 



x)l 



j=-oo 



+ ^T{2-'-)F^:l{k,^ = Q■a) 



(A12) 



Now, perform the derivative with respect to k at k = 0. Define I^gp — gp- 

expanding the argument of the F- functions in e(= 4 — d), employing the representations 



Eqs. f[T2|) and ([T3|) to resolving the remaining summation and neglecting higher order terms 
in e, we get to 



4sk^) = - 



1 

4^ 



+ 



i-y 



x{l — x)a^ 



(1-e) [ dx[x{l-x)] ^ [ dyy{l-yf2\y{l-y)T^ + ya ^ 
Jo Jo '- 

+ ^ dx[x{l — x)]~'2 / dyy{l — y)^ 
Jo Jo 



X fi.Jyr, . y{l - y)T^ + a ^y 



i-y 

x(l — x) 



-2eFi{a) 



(A13) 



This is the explicit form that can be reduced further by noticing that 0(e) terms in this 
expression can be computed at e = 0. We can simplify the final steps when we write last 
equation in terms of the following parametric integrals 



—2 f dxdyy 



In 



yi^ - y)r'r' + y + 



1 

2' 



i7(^)(r) = 2/ dxdyy x 
/i (y^, \l 2/(1 - yW + ^"^2/ + ^^^(1=^ 



(Al4a) 



(A14b) 



with r = a. Inserting these definitions in the remainder of last integral we finally obtain 

(A15) 



where W^'^\r) = G'^'^\a) + H^'^^a) — 4Fq ^(cr). Let us compute the three- loop contribution 
for the two-point function following the same reasoning. We can write that integral in the 
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form 



J=-CO 



5i 



(A16) 



[{q-kf + a\j + Ty + l\ 

Making use of the result for the one-loop bubble and neglecting higher order corrections in 
e we obtain the following intermediary result 



j=-oo 

dx 1 2 



[(g - A;)2 + a%j + r)2 + 1] 



+ eF^j\q,(r) I - + j V + 1]^ .(A17) 



(1-6) X 



2 I ■2_2 I 



Iq [x(1 — x)g2 + ^27-2 _|_ 2 

It is not difficult to prove that we can rewrite last equation as 

^ 00 
4^\k,z = 0;a) = -\a{l + e) 



dx 



J = -CO 



+ sM^l{k,a) + 0{e') 



d'^-^q 



[{q - kf + + r)2 + 1] 



[^2 + or2j2 + 

Define the subintegral 

00 

^{k,cr,x) = ^ 



(A18) 



(A19) 



- [il - k? + a\3 + r)2 + 1] [q^ + a2j2 + 

We utilize another Feynman parameter in last integral in order to compute the momentum 
integral employing Eq. (fTOj) . The result in terms of parametric integrals reads 



i5(A;,cr, x) 



g,,ir(^)r(-i + | ^ ri 
2r(e) 



dy{l - y) 



^ y(l - y)[e + aV^] + (72(j + yrf + y + 



j=-oo 



1-y 

x{l — x) 



(A20) 



Once again, use the identity y/T:T{^H^)Sd~i = V{^)Sd and expand in e. Replacing it back 
into I^^\k,a) and absorbing the factor Sd, we find 



4^\k,t = 0,a) = ^\a 



ir_r(2 -t)r(-i + |) 



X 



00 ^ 



_, 1_3£ 

2 2 



j=-oo 



+ eFi^Uk,i = 0;a)\. 
2 ' J 



(A21) 
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Now performing the derivative and defining I'^gp = ^^^-4^^-^|fc2=o results in the expression 



/'M (a) - - 



1 3e 



a 



r(2-i)rH + ¥) 

2v^r(e) ^■ + ^1^2 2 



dx 







X 



X J 



j=-oo 



(A22) 



Performing the summation we can rewrite this integral as 



r(2-f)r(-i + | 



2r(i + f)r(e 



1 3e 

2 ~ T'* Jo Jo ~ 



X 



2 L a; 1 — a; o"^ 



3e 
2 



i+^lyr, Jy(l-y)r2 + y(7 2 + — — 

2^ 2 V y x{l — 



(A23) 



Expanding the F-functions and setting e = in the subscripts of the functions which remain 
just as we proceeded in the calculation of the two- loop contribution we obtain the result 



(A24) 



The integral Iasp{<^) can be solved along similar steps. Realizing that, it can be written 
at the symmetry point in terms of the four-point one-loop subdiagram as 

j=-oo 

we can solve the internal bubble, and obtain 



(A25) 



1 pi ^ p 



j=-oo 

1 



(A26) 



[(g2 + cr2r2)x(l-x) + l]2 ' 2' t- 

Before proceeding, let us prove that last term is convergent, does not contribute to the 
ultraviolet divergences of I^^^picr) and therefore can be neglected in the consideration of its 
singularities (dimensional poles in e). Explicitly, it is given by Eg. (1A6I) 



spoo 



d'^-^q- 



(A27) 
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Since this is a massive integral, the potential singularities come from the region of high 
momentum. From simple power counting together with our previous dicussions, this integral 
will be divergent if, in the limit q — )■ oo, the F^\q,j; cr) behavior is or proportional to 

2 

a positive power of q. Thus, it suffices to prove that this object is proportional to a negative 
power of q, which we shall show next. From the definitions given in the main text Eq. (ll7p 
and the representation in terms of sum involving the product of cosine and Bessel function, 
we can write it as 

„1 oo c 

Fi'^\q, i; a) = Aa^' dxY^ cos\2nm(T + ix)] ( =) ' 

(^27rma~Va;(l - x){q^ + aH"^) + ij. (A28) 



X 



In order to attain maximal simplicity, take z = in the above expression, since it is obvious 
that F^\q,i; a) < F^\q,i = 0;a). Note that the integrand is symmetric around x = ^ 
and we can write 

(g, i = 0; a) = 8a-' dx^ cos(27rmr) ™ 



m=i ^ ^^x{l-x)q^ + l' 

X K^(2Txma-^^x{l-x)q^ + l). (A29) 

In the limit g — t- oo, choose a small real parameter A << 1 with the property \q^ — )■ oo. 
The idea is to split the integration limits into two pieces: in the first one we use the Bessel 
function and in the second piece we replace its asymptotic form for large values of the 
argument, namely 



limF^\q,i = 0; a) = 8a^^ / dxcos{2n'mr) ( 

m=l ''^ 



nm 



^(7 xq^ + 1 



,1 

~"" ^ 2 



X K^(^27ima ^\/ x(f + 1 j + ^ dxcos{27imT)(^ I x 



~(7^^\/x{l — x)q^ 



TT 



Anma ^ x{\ — x)q^ -\- 1 



exp[ —2'n-ma ^y^x{l — x)q'^ + 1 



(A30) 



The second term can be neglected in this limit. Performing the change of variables y = 
1 + xq"^, we can rewrite last expression taking into account these observations in the form 

(t) / \ ^ V ■\ , , , _t 

limFl {q,i = 0;a) = — — > cos(27rmr)(7rm) 2 / dyy 4 
q^oo 2 q^ , Jl 

x K^{2Tima'^y). (A31) 
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Using the identity [22 1 

/oo 

we are able to prove that 



u—fj.yO') 1 



limF^\q, i = 0;a) 

q—^oo 2 



^cr 2 



cos(27rmr)(7rm) 2 ^ir|_i(27rmcr 



(A32) 



(A33) 



m=l 



which is clearly regular in e completing our task in proving that in the ultraviolet region 
this object and therefore the desired integral involving it are both finite. Then it is safe to 
neglect that term in the computation of I^Jp{a). 

When we neglect the correction from Eq. flA26l) as just discussed, we can proceed from 
the following expression 



„i _ °° r 

-(1--) J dx[x{l-x)]^ J2 J 



j=-oo 



[g2 + a2(j + r)2 + l]^ 



We introduce an additional Feynman parameter y to rewrite last equation as 



(A34) 



lSp{a) = a 



r(2 + f) 
r(f)e 



„1 oo „i 

dx[x{l-x)]^ V] / dyy{l-y)^2-^ 



J=-oo 



d'^-\ 



2+f 



(A35) 



g2 + o-2(j + yr)2 + y{l - y)a'^T^ + y + 
Utilizing Eq.f llOp . we resolve the momentum integral and express the result in the form 



r(i + f)s,„ir(^^) 



^r^:y^ ' ' ' (l-^) I dx[x{l-x)]^ I dyyil-yy^-' 



j=-oo 



X) 



(A36) 



Performing the summation using Eq. f|T2l) . employing the identity \/TTT{^^)Sd-i = T{^)Sd, 
expanding the argument of the F-function in e and absorbing the factor 5*^, we obtain the 
following result 



-2e 



F(2 



2F( 



' dx[x{l — x)] 2 / dyy{l — y)2 ^ 
'0 Jq 



Tie){yil-yy + ya 



-2 , (l-y)^" 



x(l — x) 



+f^+AyT, Jy{l -y)T'^ + ya- 



x(l — x) 



(A37) 
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Note that the integral over y possesses a pole m.y = 1. A standard procedure is to compute 



16|, which facilitates the computation 



the bracket which multiplies the integral ai y = 1 
and retains the pole contribution which we are interested. Expanding in e and neglecting 
0(e°) terms, we finally obtain 

1 



2e2 2 

These integrals and all the nomenclature defined here are utilized in Sees. II and III. 



(A38) 



Appendix B: Massless integrals in dimensional regularization 

The higher-loop massless integrals are even simpler to evaluate than those occurring in 
the massive setting discussed above, with the difference that they need to be calculated at 
nonvanishing external momenta. The massless counterparts of the integrals discussed in the 
previous Appendix are given by 

d'^-^qid'^-^q2 



1 



(g? + ^^(ji + r)2)(g2 + a2(j2 + r)2) 



[(gi + g2 + A;)2 + a2(ji+j2 + 2 + r)2] ' ^^^^ 



f d qid q2 



Jl j2 = -00 

1 

{ql + ^'(J2 + r)2) [(gi - + k^)' + a^{j, - J2 + ^3 + r)'] 
1 

[{P - qi)^ + c^{p - Ji + t) 



00 „ 

l5{k,i,a)=a^ ^ / 



d'^-^qid'^-^q2d'^-^q3 



{q! + cy\ji + ry)iql + a\j2 + ry) 
1 



{ql + ^'Us + rf + 1) [(gi + q2 + kf + a^j, +32 + ^ + rf] 

""[{qi + qz + kf + a^n+n + i^rYX ^^^^ 

Firstly, let us discuss Is at arbitrary external momentum, which shall be important in the 
evaluation of critical exponents using minimal subtraction. As we are going to compute the 
exponents in normalization conditions as well, we shall discuss its derivative with respect 
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to squared nonvanishing external momenta for that sake. All integrals will be computed 
at zero external quasimomenta from now on. Since many developments discussed in the 
previous Appendix will have their exact analogy here, we will be more economical in the 
steps to compute the integrals. First, write in terms of the one- loop subdiagram as 

4-\k,z = 0;a) = a V / d'-\- ff,^^,^'^',^^ (B4) 

J ^ [{q - ky + + ry] 

Replacing the value of the subdiagram computed in Sec. IV through the use of Eqs. fH5]) - flSU]) 
defined for the massless case, we get 

i^j\k,i = Q-a) = r'c/xb(i-x)i-t r . 

+f r(2 - f )^ / t(,_,)2^:;;/,.)2^,] i^^^ (g, I = 0; a)) . (B5) 
For the sake of minimal subtraction, define the parametric integral 

L3(P + aV^) = /' dyylnlvil - y){k^ + a^r% 
Jo 

Other useful definitions are the massless functions 

p(r)., . s_ 1 f ,d-i Ft\q + k,j 

Let us compute i^{k,a)^ which is analogous to the expression for the massive case and 
much simpler. We introduce a Feynman parameter y in order to express the two denomina- 
tors as a single one, solve the momentum integral using Eq.f llOp and expand in e afterward. 
This set of steps are identical to those which led to Eq.f lAllj ) and the reader can check that 
it results in the following expression 



(B6) 



H{k,(T) = dy[l-y)'^ 

oo 1 



-1 



- y){k' + aV^) + a\3 + yrf] . (B9) 

=— oo 



j=-oo 
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The summation can be performed as before and we find 



i3{k,a) 



dy{l - y)-^-'V^ a-'T{e - 1) y{l - y)[k 



2r(|) ,0 



(BIO) 



Utilize the identity ^/7^^{^)Sd-l = T{^)Sd, absorb in the redefinition of the couphng 
constant just as explained for the one-loop four-point graph in the main text and expand 
in e. Integrating over the Feynman parameter x and substituting into the expression of 
I^^\k, i — 0,a), we obtain 

) {k,t = Q;a) = -({l- 11 /' dyil - yf^-^ [r(e - 1) 

3 V > ' ^ 2' 2r(|)r(2-e) Jo y 

l-e 



X (y{l - y)[e + aV]) " V^-^7-i+e(2/^, a'' ^ y{l - yW a^^]) 



(Bll) 



Next, use the property r{x + 1) = xT^x) in order to get a result purely in terms of e when 
higher order terms are neglected. We then have 

l!f\k, ^ = 0; a) = 1 (-1(1 + e) dyil - y)i-^ [(y{l - y)[k' 

+ '-F^Jl{k,t = 0,a)). (B12) 

We introduce two new parametric functions which will be useful to our future manipulations, 
whose expressions are given by 

Fi^\k, i = 0; a) = a^-^- dy{l - yT^-^f.^^yT, V^(r^^)(FT^),(B13a) 
Fi'\k,i = 0;a) = dyy{l - y f^f^.^yr, a-'^y{l - y){k' + aH^)). (B13b) 

In terms of these functions, we can express the solution for the integral l3{k, i = 0; cr) in a 
form appropriate to minimal subtraction, namely 



8e 



AeFi^l{k,z = 0;a)). 



(B14) 



If we prefer to employ the normalization condition scheme, we need the derivative of this 
integral with respect to k^ computed at the symmetry point k^ — 1. Henceforth we denote 
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the argument of a function {k,i — 0,a) at — 1 by (a). Using the recursion relation for 
Bessel functions 

dlz-'^K^iPz) 



dz 



-/Sz-'^-'Ko^+iiPz), 



(B15) 



and employing the explicit representation of /a(a, b) in terms of the product involving the 
Bessel function and cosine, we learn that 

iSU'^) = -^(1 + T - 2^^^(1 + ^V) + 26Ff)(<7) - ieF'f\a)). (B16) 

Let us conclude our analysis of the two-point function by calculating the three-loop 
integral l'[\k,i — Q;a). First, write it as 

]=-oo 



A;)2 + a2(j+r)2]' 



(B17) 



When we use the explicit expression for the four-point one-loop subdiagram along with 
the value of the integral over the parameter originally appearing there, expanding the F 
functions in e and neglecting higher order terms, we obtain 



4^)(A;,z = 0;a) = ^(a(l + 6)Er=-oo/ 



[(9-fc)2+<T2(i+T)2][g2+a2j2]. 



);<7)). 



(B18) 



Using completely similar steps in the calculation of llf\k, i = 0;(t) with minor modifications, 
we can proceed henceforth quite analogously Define the functions 



F^^\k, t = 0;a) = a'-'^ dy{l - yf-' f_.,.AyT, a-'^y{l - y){k^ + aH^)), (B19a) 

2 2 2 

F^Z\K i = 0; <j) = a-^^ dyy{l - yYh+s.{yT, a'Wvi^ - vW + 'j't^))- (B19b) 



The solution for useful in utilizing minimal subtraction scheme can be written as 



lt\k,a) 



6e2 



({k^ + C7V)[1 + J - 3eL3(A;2 + aV^)] - 3eF^Z\k, ^ = 0; a 

2 



(B20) 



Instead, if we wish to employ normalization conditions, we have to compute the derivative 
in relation to k^ at the symmetry point k'^ — 1 and we get to 

4p(^) = -^(l + 2^ - ^'Hh^ + '^'^') - + ^(^)])- (B21) 
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It is worthy to mention at this point that the finite size corrections for l'^sp(cr) and 
I'^sp{ct) become simpler when computed at e = 0. In fact, defining the massless quantity 
Wq^^((j) = + cr^r^) — Flf\a) + 2FQ'^\a), those integrals can be written simply as 



4p(^) = -^(1 + T - ^^Wo(a)), (B22a) 
r'ir) 



((7) = --^(l + 2e-3eW^o((7)). (B22b) 



The massless counterpart of the integral which contributes to the four-point function at 
two-loops can be written in terms of I2 and reads: 

Once again we shall restrict ourselves to the simplest expressions for this integrals, namely 
those with zero mode (p = 0) characterizing the external quasimomentum associated to 
the finite size direction, perpendicular to the plate surfaces. At the symmetry point, all 
primitively divergent vertex parts depend on only one external momenta scale. For this 
reason we are going to list the last integral at a certain external momentum P and display 
the results in the most convenient form either using minimal subtraction or normalization 
conditions renormalization schemes. 

We solve the internal bubble, i.e., we compute the integral over the momenta. Then, use a 
Feynman parameter x to melt the integer powers of the propagators in a single denominator. 
We are left with an integral over the Feynman parameter x as follows 

&\ki,0;a)^[-(l + -)a dx Y ^ ^7 ^ 

^ ^ ^ ^e^ 2' Jo J [g2-2a;Pg + a;P2 + cj2(j+r)2]2 



j=-oo 



Xjr—Tr4—-^ + l^fiP^ 0, a)Sai (B24) 



where 



rir)(pj, n^\^^\^ f d'^ ^qF^{q + k3,j,(T) 

" ^ ' ' ^ S.p^J [{q-Pr + a'{j + rm' + <^'U + rrY 



(B25) 



j=-oo 

We shall be concerned with the behavior of Fa{q,j,a) in the limit g — >■ 00. In that 
limit, if this object behaves as an inverse power of q all is well, since last integral will have 
no dimensional poles in e using simple power counting arguments. Therefore, it can be 
neglected in our consideration of poles from Ii{k,i — 0, cr). Let us turn our attention to this 
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issue. From its explicit form Fa{q,j,a) < Fa{q,j = 0,0") and we just have to discuss the 
latter in the appropriate limit. 

From our discussion in the main text Sec. IV, we already know that 

( ) e e °° el 

FJ {k, z = 0; 0") = 4A;~5(Tt ^ cos(27mr)(7m)5 dx[x{l — x)\~^ 

2 n=l 

xKL{2Tinka-\x{l-x)]^. (B26) 



The integrand is symetric around x = |. Then we can multiply the integral by two and 
use the integration limits at (0, |). Next, perform the change of variables a;(l — x) = . 
Thus, we can write the integral in the form 

f dx[x{l-x)]-'iK^{2'Knka-^[x{l-x)]^) = 2^ f dyy^'^l-y^)'^ K^^nnka-^y). (B27) 
Jo ^ Jo ^ 

After some manipulation with special functions using Ref. 22|, last integral can be put in 
the form 

/ dx[x{l — x)]~^ K±(27mka^^[x{l — x)]^) = 2^ -cosec( — ) cosh^nnka^^) 

r''^''~\sinht ^ , . , l^ r''''''' cosht - 1 \l 

/ dt- stnh{7mka'^)[-f + ln{TTnka'^) + / dt) . (B28) 

Jo ^ ^ Jo ^ ^ 

The first integral in the right-hand side is defined as the integral hyperbolic sine and denoted 
by shi{7Tnka^^). The complete term multiplying sinh{7mka^^) is defined as the integral hy- 
perbolic cosine whose symbol is chiijmka^^) . Using the identity coseci^Y) = ;^ + f| + 0(e^), 
the above expression becomes regular in the limit e — )■ 0. Now, the finite size contribution 
is 0{e^) (regular) and we can set e = inside its expression in order to obtain 

oo 



Fi^\k,i = 0;a) = \cosh(nnka-^)sM(7rnka-') 

n=l 

-sinh{-nnka-^)chi{-nnka-^)] + 0(e). (B29) 



Consider the terms inside the brackets and define z = imka ^. For example, from our 
definitions, we can write 

coshzshiz — sinhzchiz = —{shiz — chiz) H — —{shiz + chiz). (B30) 

We can work out these equations further in order to reduce them in terms of the incomplete 
Gamma function defined by r(0, z) = dt^, such that 

coshzshiz — sinhzchiz = -—T{0,z) H — ;^r(0, —z). (B31) 
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Now the limit A; — )■ oo is the same as z — )■ oo. From Ref. 
is given by limT{0, z) = Therefore, we find 



22| . the asymptotic value for real z 



limicoshzshiz — sinhzchiz) = —. (B32) 

z-> z 



Thus, we conclude that 

limFr{k,i = 0;cr = —3^^^, B33 



n=l 

showing that this behavior guarantees that the integral G^[\k,i = 0, a) is regular and can 
be neglected in the determination of the coefficients of the pole terms contained in this 
diagram. 

In the following we outline the method of computation of I^\k,i = 0,a). It does not 
bring essential new features in comparison with the massive case, being actually simpler 
its evaluation. Consider the first term in the last expression l[^\k,i = 0,a). Using an 
additional Feynman parameter y to fold the remaining terms into a single denominator in 
order to compute the integral over the momenta q, the resulting parametric integral over y 
has a pole in y = 1. In analogy to the massive case in Appendix A, we keep the prefactor 
(depending only on y) in that integral, and set y = 1 in the overall term which depends on 
(x, y, P, a). We then find 

li^\h,0;cr) = ^(1 - I - ^^iP) + eFW(^,0)), (B34) 

which at the symmetry point = = 1 takes a simpler form and simplifies further if the 
finite size correction is computed at e = 0, namely 

4U^) = ^(l + | + e^o^^H^))- (B35) 

We are going to use these results in Sections IV and V in the massless computation of the 
critical indices. Notice that the results are such that F!i'\a) used in the present appendix 
is the appropriate object for the massless theory, given in the text by Eq. (H9!) . 
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